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EFFECT OF FACE -SHEET STIFFNESS ON BUCKLING OF 
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SANDWICH CONSTRUCTION IN 

AXIAL COMPRESSION 

By Robert E. Fulton 
Langley Research Center 

SUMMARY 


A study is made of the effect of the flexural stiffness of the face sheets 
on the buckling of elastic curved plates and cylindrical shells of sandwich con- 
struction subjected to axial compression. The study shows that when the core is 
very weak in shear, the flexural stiffness of the face sheets can have an impor- 
tant effect on the buckling load. Simple formulas are developed which give good 
approximations to the buckling load for most practical ranges of the parameters 
where the effect of face-sheet stiffness is important. Results of these formu- 
las are compared with exact results for an infinitely long curved plate over a 
large range of sandwich- shell parameters. The results are based on a linear 
buckling theory for sandwich shells of the Donnell type (NACA Report 479), which 
takes into account the asymmetry of the sandwich faces. 


INTRODUCTION 


Over the past few years a number of studies have been made of the buckling 
characteristics of curved plates and cylindrical shells of sandwich construc- 
tion subjected to axial compression. In most of these studies the contribution 
of the flexural stiffness of the face sheets to the buckling characteristics has 
been neglected. Although the assumption that face-sheet contribution should be 
neglected is reasonable in certain ranges of the significant parameters, it may 
not be reasonable in the range where large shearing deformations occur in the 
core at buckling. Buckling in this latter range of parameters results in small 
axial buckle wavelengths; consequently, the faces undergo large bending deforma- 
tions. The strain energy associated with this bending may be significant, and 
the effects of face-sheet stiffness should be considered in determining the 
buckling load of the shell. 

The purpose of this paper is to investigate the contribution of the flex- 
ural stiffness of the face sheets to the buckling load of a curved plate or 
cylindrical shell of sandwich construction in axial compression. The investi- 
gation shows that a simple expression for the buckling load including this 



face-sheet contribution can be developed which holds for most practical ranges 
of the parameters where the face-sheet contribution is important. The present 
study is an extension of the work carried out by Stein and Mayers (ref. 1) where 
the face-sheet stiffness was neglected. 


SYMBOLS 


a,b 

B 


coordinates of edges of shell in x- and y-directions, respectively 

Et 


extensional stiffness of face sheet. 


thickness of core 


flexural stiffness of face sheet, 


1 - 


Et^ 


12(l - n 2 ) 


d = 


d = 

E 

F 

Gc 

h 

ka 


k b 


cB-j_B2 

G c (Bi + Ba) 

h 2 B-LB 2 

Bi + B2 

Young* s modulus for face sheet 
Airy stress function 
shear modulus for core 


distance between middle surfaces of face sheets 


buckling-load coefficient for cylinder, 


N*a 2 ( B i + B g ) 
rt 2 h 2 B 1 Bg 


buckling-load coefficient for curved plate, 


M*b 2 (B 1 + B 2 ) 
:rt^h^B]_B2 


m,n integers 

M x ,My sum of bending moments in face sheets 

N cr critical axial stress resultant 
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buckling increments in normal and shear stress resultants 




externally applied axial stress resultant 
externally applied normal and shear stress resultants 


radius of curvature of shell 


face-sheet-stiffness parameter. 


l B i + b 2 )(Di + d 2 ; 
h 2 B 1 B 2 


thickness of face sheet 

buckling displacements of shell in x-, y-, and z-directions, 
respectively (fig. l) 

constant (eq. (6)) 




coordinates of middle surface of shell (fig. 1) 


I o a^B-i + Bp) 2 (l - n 2 
curvature parameter for cylinder [Z a = — — * * 

\ r^h^B^B2 


curvature parameter for plate IZ^ 


+ B 2 ) 2 (l - 4 2 ) 




buckling rotation in x-direction, ^( u l “ u 2) 
buckling rotation in y-direction, ^( v l “ v 2 ) 


length -width ratio, 


Poisson’s ratio 


coupled rotation variable, a, x + p,y 


sandwich-core parameter for cylinder. 


3t^cBj_B2 
i 2 G c (b ± + B 2 ) 
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sandwich-core parameter for curved plate , 


Tt^cBpBp 




b 2 Gc^B-]_ + B2^ 


V 2 


two-dimensional Laplacian operator 


Subscripts: 

i integer 1 or 2 

1,2 refer to upper and lower face sheets, respectively 

x,y after commas, indicate partial differentiation with respect to 

axial and circumferential coordinates, respectively 


GENERAL BUCKLING EQUATIONS 


Assumptions in the Sandwich Theory 

The buckling equations to be solved were obtained from the nonlinear equa- 
tions originally developed in references 2 and 3- The variables in the nonlin- 
ear equations were separated into prebuckling and buckling effects in a manner 
similar to that in reference A. In the present study boundary conditions on 
the prebuckling state have been relaxed so that the prebuckling displacements 
are either constant or linear. 

The concept of a sandwich is retained, in that the core undergoes only 
transverse shear deformations so that a line through the undeformed core remains 
straight when the core is deformed but does not necessarily remain perpendicular 
to the neutral surface of the shell. It is assumed that the total thickness of 
the shell element is small compared with the radius of curvature. The face 
sheets are assumed to be elastic, isotropic, and homogeneous and to follow clas- 
sical shell theory - that is, to satisfy the Kirchhoff-Love condition. The core 
is assumed to be elastic, isotropic, and homogeneous, to carry no inplane loads, 
and to have no deformation in the direction normal to the neutral surface of the 
shell. Poisson’s ratio is taken to be the same for the two face sheets and the 
core; however, the thicknesses and Young's modulus may be different. 


Governing Differential Equations 

If the above assumptions are employed, the three differential equations 
governing the buckling of a cylindrical shell of a nonsymmetrical sandwich sec- 
tion (fig. 1) are as follows: 

+ B^l - (la) 


k 



(l - av 2 )*? = v 2 v 


(lb) 


(l>l + D 2 )v^w - N*w ;xx - 2N^ r w, xy - Njv,yy - i F,^ + avScp = 0 (1c) 


where 


T X ^ >y 


d = 


cB^Bg 

G C (Bl + Ba) 


a = 


h^B]_B2 
B^ + 


a . i(u L - us) 6 - Wn. - v 2 ) 




D-i = 


Eiti 


.3 


12 


t 1 - n 2 ) Bl " t 1 - 1* 2 ) 


E^i 


(i = 1 or 2) 




( 2 ) 


The quantities u-j; , Vi , and w (wi = w) are the middle -surface incremental 
buckling displacements of the ith (upper or lower) face sheet of the sandwich 
in the longitudinal (x), circumferential (y) , and inward radial (z) directions, 
respectively. Young’s modulus for the face sheets is denoted by E^, Poisson’s 
ratio by \±, and the shear modulus of the core by G c . The thickness of the 
ith face sheet is given by ti, the thickness of the core by c, the distance 
between face-sheet centers by h, and the radius of curvature of the shell by r. 
The subscripts x and y after a comma indicate differentiation with respect 
to the axial and circumferential coordinates, respectively. The Laplacian oper- 
ator is given by V 2 , and F is the Airy stress function such that the buckling 
increments in the normal stress resultants N x and N y and in the shear stress 

resultant Nxy are given by 


N X = F >yy 

% = f >xx y 

^xy = _F >xy 


(3) 


The inplane forces which exist in the shell just prior to buckling are denoted 
by N^, (assumed to be positive in tension) , and Nxy. 


The buckling behavior of the sandwich shell can be determined from a solu- 
tion to equations (l) with the appropriate boundary conditions in terms of the 
three unknowns, w, F, and cp. The usual shell theory, which includes the 
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effect of shearing deformations, requires that five conditions he satisfied on 
the boundary. For equations (l), however, six conditions are necessary. The 
extra condition results from both including the flexural stiffnesses of the face 
sheets in resisting the total moment and retaining the Kirchhoff-Love assumption 
for the individual faces. Thus, although in classical shell theory only one 
boundary condition is required for the total moment normal to an edge, the pres- 
ent study requires two separate conditions, which separate the -total moment into 
a moment in the faces and a moment couple resulting from inplane forces in the 
faces. 

In equation (lb) the variable cp results from coupling, through differen- 
tiation and addition, the two equations which define a -and 0, the rotations 
in the longitudinal and circumferential directions, respectively. (See ref. 2.) 
A complete set of boundary conditions for equations (l) should therefore be 
specified in terms of a and 0, not just cp. In the present study, however, 
it is convenient to use equation (lb) because the boundary conditions are such 
that cp, rather than a and 0, is a natural variable. 

BUCKLING OF A CURVED PLATE 


Plate of Finite Length 

The problem of interest is to determine the buckling load of a rectangular, 
simply supported, cylindrical plate subjected to a normal compressive force N* 

parallel to its directrix along the edges x = 0 and x = a (fig. 2). For 

* -x- -x- * 

this loading condition, N x = -N and N xy = Ny = 0. 

The six boundary conditions which govern buckling and which are defined 
herein as simply supported are given as follows along x = 0 and x = a: 

(l) Displacement normal to the surface of the plate vanishes; that is. 


w = 0 (4a) 

(2) Moment couple normal to the edge due to differential inplane forces in 
the face sheets vanishes; that is. 


a 'x + y = 0 


(4b) 


(3) The sum of the moment in each of the individual face sheets vanishes; 
that is. 


My = 


= -(Dl + P 2 )(w, xx + = 0 


(4c) 
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(4,5) Displacements parallel to each edge are prevented; that is. 


v 


B l v l + B 2 v 2 = Q 

Bi + B2 


(4d) 


P = 0 


(4e) 


(6) Motion normal to each edge in the plane of the sheet occurs freely; 
that is. 


F ’yy ‘ 0 


(Uf) 


A comparable set of conditions along y = 0 and y = b are: 


w = 0 


P>y + = 0 


Mt 


- -(D! + B 2)^ w ^yy + H^xx) “ 0 


u = B l u l + B 2 u 2 = q 
B}_ + 


(5a) 

(5b) 

(5c) 

(5d) 


a = 0 ( 5e ) 

= 0 (5f) 


It should be noted that equations (4b), (4e), (5b), and (5e) imply that cp 
vanishes on the boundary. The quantities M x and My are the sum of the 

moments in the two face sheets and u and v are the neutral surface displace- 
ments in the x- and y-directions, respectively. 

The solution to equations (l) is 


w = w 0 sin 322 sin 22£ (6) 

a b 
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F 


- (Bl ♦ Ba)(l - ^ 




” 2a2 s + s 


sin ^ sin 
v 2 a b 


(7) 


o/ m 2 n 2 , 

w o rt ~p + “p I 

V a b / • mitx . nny 

i ' — sm ==i=. sm — — 

1 + ” 2d (S + S) 


cp = 


( 8 ) 


where w Q is a constant and m and n are integers. Equations (6), (7)) 
and (8) satisfy explicitly all the applicable boundary conditions (eqs. (4) 
and (5)j except those given by equations (4d) and (5*3-) • It can be seen that the 
conditions given by equations (4d) and (5d) are also satisfied by noting that u 
and v are related to w and F by (ref. 2) 


■\ 


u 


F 'yy " 


xx 


’X 


(l - u 2 )(b! + B 2 ) 


F,xx - HF, 


y (l - h 2 )(bi + B 2 ) 


yy + v \ 

r / 


-2F, 


>y v 'x 


xy 


(1 - + B 2 ) 

Substituting equations ( 6 ), ( 7 ), and (8) into equation (lc) yields 

^ £ ♦ (* + d 2 y * ^X 1 - u f k 

a 2 ' ' r 2 (m 2 + A 2 jt 2 ) 2 


(9) 


cbA(m 2 + A^n 2 ) 


i4 n + S!^(m 2 + X 2 n 2 ) 
_ a 2 


w 0 = 0 


( 10 ) 
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The trivial solution corresponding to the unbuckled state occurs for w Q = 0, 
and the buckling load is obtained by setting the braced term equal to zero. 
The equation for the buckling-load coefficient is therefore 


where 


SI + 


2.2 1 


k b = 


An 


m^ 


m2 


1 + 


■\2 2' 





h ♦ 

2£ 

1 + 

m 2 1 

l * 2 ) 

m 2 



*\ 2 2 
A n 




k-K = 




N*b 2 ( B 1 + B 2 ) 
rt 2 h 2 B 1 B 2 

2 = ^ 4 (Bi + B 2 f(l - n2) 
r2h 2 B 1 B 2 


ilr = *£§ - * cB 1 B 2 

13 v2 


b c + %>) 


S = 


+ D 2 ^B^ + B 2 ^Dj^ + P 2 ^ 


h 2 B 1 B 2 


A = 2. 


( 11 ) 


(12) 


If the face-sheet flexural stiffnesses Dq and D 2 are neglected as in 
reference 1, equation (ll) with S = 0 becomes 



( 13 ) 
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Equation (13) is identical with an equation developed in reference 1 except that 
the core-stiffness and curvature parameters used herein are more general and 
take into account the asymmetry of the sandwich section. 

It should he noted that equation (ll) is quite general for all ranges of 
the parameters and is applicable for both curved plates and cylinders. If the 
core thickness is zero (c - 0), equation (ll) reduces to a Donnell-type (ref. 5) 
result for a bilayered shell. With c = 0 and the face sheets symmetrical, it 
reduces to the equation given in reference 6 for an isotropic shell. 

The buckling load for a shell of finite length can be obtained by minimizing 
equation (ll) with respect to m and n. For an infinitely long shell the mini- 
mization should be carried out with respect to n and to the buckle -wavelength 

ratio -2-. 


Plate of Infinite Length 

Equation (ll) has been minimized with respect to n and — and has been 

An 

applied to an infinitely long curved plate. The results are given by the solid 
lines in figure 3- It should be noted that the results for S = 0 shown in 
figure 3 agree with those in reference 1 and that the results in figure 3( a ) for 
S = l/3 and c = 0 ^tq = tg = hj agree with those for an isotropic homogeneous 

curved plate given in reference 6. 

Stein and Mayers (ref. l) give the following formulas for the critical load 
of an infinitely long curved plate obtained from an equation of the form of 
equation ( 13 ): 


z h k \P- - tb 

(1) For 

* (1 + tb) 


k b - it + 

(i + *o 2 x 1 * 11 


(14) 


(2) For W 1 - fa < Zj, < £ - % 
(i + tb f * 2 ' 


tb 




it c 


z b , 
*2 +* 


V 1 - tb\ 


( 15 ) 
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(3) For 


Zb > \/ X - *b 

* 2 = +b ' 


- £ (16) 


Although figure 3 shows plots of the critical load for all ranges of the 
parameters, it is of interest to develop simple formulas which apply for the 
more important ranges. In particular, the range where S is less than 0-01 is 
of practical importance. For a symmetrical sandwich this range corresponds to 
the case where t/h is less than about l/6. As is seen from figure 3, it is 
reasonable to consider equations (l^-) and ( 15 ) to be applicable to sandwiches 
where the face -sheet thickness is small with respect to the depth of the sand- 
wich. On the other hand, equation (l6) is not adequate even for the range 
S < 0.01; that is, the face -sheet-flexural-stiffness contribution is quite sig- 
nificant even when the thickness is small. 

It is well known that equation (l6) is used for S = 0 in the parameter 
range of 


5b> 

*2 " 



+b 


and corresponds to large core shearing deformations and thereby to an infinite 

wavelength ratio — . In order to develop an equation to replace equation (l6) 
An 

for a small value of S > 0, the contribution of the face sheets to the critical 

load of the plate in the range of large ~ is investigated. 

An 

In the study of this contribution for the infinitely long plate, n is 
set equal to 1 in equation (ll) and the equation is rewritten in the following 
form: 
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In a similar fashion, equation (13) becomes 


kb 



(18) 


It is clear from equation (18) that as the wavelength ratio — = — approaches 

A a/m 

infinity, equation (l8) becomes equation (l6). In equation (17) , however, as 

=2 becomes large the contribution of the first term and consequently the buckling 
A 

load become very large. This inconsistency for large values of S between 

A 

equations (17) and (18) reflects the effect of the flexural stiffness of the 
face sheets on the critical load. 


In order to determine the proper contribution of the face sheets to the 

buckling of the shell, equation (17) should be minimized with respect to S. 

A 

Unfortunately, such a procedure does not yield any simple results. Attention is 
therefore focused on the range of parameters where both » — and 

2 W +b 

\ » 1 so that equation (17) can be written as 


k b 




*b 


Minimizing k-^ with respect to 


m 

A 


leads to 



- 2 



= 0 


/ m\ 4 = Zb! 

W "As 
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The critical buckling-load coefficient becomes therefore 



For a symmetrical sandwich the critical force is 


( 19 ) 



The usual result, which is obtained by neglecting the faces, is 


( 20 ) 


G c h^ 

N - “-I" 

The results from equation (19) are shown in figure 3 as dashed lines when 
they differ from the exact results obtained by minimizing equation (ll). The 
plots indicate that equation (19) is quite accurate for small values of S in 
the range of interest, namely 


&b ^ v 1 - *b 

* 2 = *b 


Equation (19) is also seen to apply for large values of S if is also large 

and is in fact the asymptote for the results obtained by an exact minimization 
procedure. 

It is interesting to note that for a symmetrical sandwich the critical 
force N cr given by equation (20) is made up of the sum of three classical 
buckling forces for three different elements. These forces are - 

The classical sandwich force obtained by neglecting the face-sheet 

flexural rigidity and 

The classical buckling force for each of the two face sheets 
considered separately 

When the face sheets are not alike, the result is somewhat more complicated. 
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BUCKLING OF A CYLINDER 


Equation (ll) can also be used to determine the critical load of a closed, 
simply supported cylinder of arbitrary length. For a cylinder, however, 
b = 2xr and n is restricted to even integers. Plots similar to those in fig- 
ure 3 could be prepared for the infinitely long cylinder; however, this work 
was not undertaken in the present study since the trends would be similar to 
those for the curved plate. The major point of practical interest is the devel- 
opment of a simple expression for cylinders which indicates the contribution of 
the face sheets in the range of large values of Z a . In order to investigate 
the cylinder it is convenient to rewrite equation (ll) as 


ka 


^■2 


= sm2 i + ZiJi: 


5a! 


im 


2 1 + 


aM 


m c 


m^ 


m 2 * (l + 1 + * msfl + 

{ m 2 / ’ a ^ m 2 y 


( 21 ) 


or, if the face-sheet stiffness is neglected (S = 0), 


where 


ka = 


5a! 


m< 


i 2 1 + 


.2 o' 

A n^ 




m2h + 1 + V 1 ^ 1 + 




m^ 


(22) 


ko = 


*a 2 (B x + Bg) 
jt^h2B]_B2 


2 _ a 4 (Bl + B 2 ) 2 (l - q 2 ) 


r 2 h 2 B 1 B 2 




= 


5t 2 cB^B2 


a^G, 


; (bi + b 2 ) 


J 


( 23 ) 


and where S and A are the same as defined in equation (12), 
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Equation (22) is identical with an equation originally developed in ref- 
erence 1, except that the core-stiffness and curvature parameters again take 
into account the asymmetry of the sandwich section. Likewise, in reference 1 
the following formulas are given for the critical load of a simply supported 
cylinder obtained from an equation of the form of equation (22): 


(l) For 


z a < 1 

it 2 1 + t a 


k a 


1 

1 + +a 


+ 



(2) For ^ — <; 

1 + \(r a rt 2 t a 


ka = 




(3) For 




(2k) 


(25) 


(26) 


In order to develop a simple expression which includes the contribution of the 
face sheets , it is useful to proceed as in the previous case for the curved 
plate and investigate the behavior of the cylinder in the parameter range 


— > U. 

jt 2 ^a 


This range corresponds to large shear deformations in the core and a 


large wavelength ratio — . When 

An 


becomes 




» 1, equation (21) 


ka = Sm 2 


Z a 2 i 

Jt^m^ 


Minimizing this equation with respect to m 2 results in the following expres- 
sion for the critical load in the region where face-sheet flexural stiffness is 
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significant 



(27) 


As might he expected, this formula is similar to equation (19), which was 
obtained for the curved plate. 


CONCLUDING REMARKS 


A study has been made of the effect of flexural stiffness of the face 
sheets on the buckling of elastic curved plates and cylindrical shells of sand- 
wich construction subjected to axial compression. The results of the study for 
curved plates are given in graphical form. The results agree at one extreme 
with classical sandwich theory and at the other extreme with isotropic shell 
theory. Any face-sheet flexural stiffness is seen to increase significantly 
the buckling strength of these sandwich shells in certain ranges of the shell 
parameters where the core is very weak in shear. In order to obtain any 
meaningful results in these ranges, the effect of face-sheet flexural stiffness 
must be included. Simple formulas have been developed in the paper for the 
buckling load of curved plates and cylinders which take into account this face- 
sheet contribution in the range of parameters of practical interest where the 
face-sheet thickness is small. 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., January 28 , 1965. 
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